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SUMMARY  PAGE 


THE  PROBLEM 


--The  Naval  Biodynamics  Laboratory  in  New  Orleans  is  involved  in  a  con¬ 
tinuing  effort  to  describe  the  kinematic  and  physiological  response  of  ana¬ 
tomical  segments  of  human  volunteers  to  acceleration  environments.  The  scope 
of  this  presentation  is  to  discuss  some  of  the  problems  that  arise  in  the 
measurement  of  the  three  dimensional  dynamic  response  of  human  volunteers  to 
relatively  short  duration  (250  millisecond)  acceleration  profiles. 

FINDINGS  _ _  -  •  - 

v  -With  pre-selection  of  accelerometers  and  careful  calibration  procedures 
six  accelerometer  configurations  can  provide  results  which  agree  well  with 
photo-derived  results  at  the  velocity  level  in  the  time  window  of  interest 
(250  milliseconds)  for  the  acceleration  profiles  investigated  at  the  NBDL. 

7 

RECOMMENDATIONS  "" 

-Accelerometer-derived  kinematic  variables .are  best  at  the  acceleration 
level  and  photo-derived  variables  are  best  at  the  displacement  level  and 
thought  should  be  given  to  combining  results  of  these  two  measuring  systems 
into  one  consistent  set  of  kinematic  variables  from  the  acceleration  to  the 
displacement  level. 

/> 
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PROBLEMS  OF  MEASUREMENT  IN  HUMAN  ANALOG  RESEARCH 
Leonard  Lustick,  Gerald  Williamson,  Marjorie  Seemann  and  James  Bartholomew 

INTRODUCTION 

The  Naval  Biodynamics  Laboratory  in  New  Orleans  is  involved  in  a  continuing 

effort  to  describe  the  kinematic  and  physiological  response  of  anatomical 

* 

segments  of  human  volunteers  to  acceleration  environments  (1,2,3).  The  scope  of  this 
presentation  is  to  discuss  some  of  the  problems  that  arise  in  the  measurement 
of  the  three  dimensional  dynamic  response  of  human  volunteers  to  relatively 
short  duration  (250  millisecond)  acceleration  profiles. 

At  the  outset  of  the  program  it  was  recognized  that  redundant  measurement 
systems  wouiu  be  required  to  measure  all  levels  of  the  dynamic  response 
(acceleration  to  displacement)  with  the  desired  accuracy  and  to  check  the  results 
of  each  measurement  system.  Initial  studies  had  Indicated  that  accelerometer 
derived  kinematic  variables  would  be  the  best  source  of  angular  and  linear 
acceleration  components,  whereas  the  angular  orientation  and  linear  displacement 
components  which  are  obtained  directly  from  photo  data  (rather  than  by  integration) 
would  be  most  accurate  for  these  latter  variables.  Further,  if  special  care  is  used  in 
the  selection  and  calibration  of  the  accelerometers  (4)  angular  and  linear  velocity 
components  derived  from  accelerometer  data  will  compare  favorably  with  those  derived 
from  photo  data  in  the  time  window  of  Interest  (250  milliseconds). 

* 

Numbers  in  parentheses  indicate  References  at  end  of  paper. 
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In  keeping  with  this,  the  dynamic  response  of  anatomical  segments  at  the  Naval 
Biodynamics  Laboratory  Is  measured  with  accelerometer  configurations, 
clnephotography  and  gyros.  The  accelerometers  and  cine photography  techniques 
are  complete  systems  allowing  the  three  dimensional  motion  of  a  rigid  body  to  be 
described.  The  gyro  configuraton  only  measures  two  components  of  angular 
velocity  and  is  used  for  monitoring  the  experimental  runs  and  for  validation 
where  applicable  (5,6). 

It  was  always  the  intention  to  finally  combine  the  measurements  from  the 
redundant  measurement  systems  in  an  "optimum'*  fashion  and  consequently  represent 
the  dynamic  response  at  all  levels  (acceleration  to  displacement)  in  one  consistent 
"best"  set  of  data.  Subsequently,  preliminary  results  of  a  method  to  accomplish 
this  optimum  combination  are  presented.  A  more  complete  description  of  the  method 
and  its  utilization  will  be  presented  in  a  future  paper. 

PHOTO  PROBLEMS 


The  target  configurations  at  the  Naval  Biodynamics  Laboratory  are  attached  to  a 

rigid  body  and  are  tracked  by  at  least  two  cameras.  The  displacement  and  orientation  of  th< 

rigid  body  are  determined  with  a  least  squares  algorithm  using  the  target  locations 

In  the  film  plane  as  well  as  the  measured  location  and  orientation  of  the  cameras  In  the 

laboratory.  At  the  NBDL,  the  coordinates  of  the  targets  in  the  film  plane  are 
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determined  by  an  automatic  digitizing  system.  The  system  in  general  is  good,  but 

editing  is  required  to  remove  macroscopic  errors  such  as  misidentification  of 

targets  and  crossing  targets.  It  has  been  found  useful  to  plot  the  target  configuration 

determined  by  the  film  coordinates  as  seen  from  each  camera  (25  frames 

per  page)  to  locate  and  to  correct  if  possible  these  macroscopic  errors.  As  an 

illustration  Figure  1  shows  25  frames  of  a  particular  run.  In  addition  to  editing, 

this  method  allows  a  quick  look  at  the  motion  as  seen  by  the  camera. 

To  process  photo  data  it  is  required  to  have  a  transformation  from  points 
in  the  object  plane  to  points  in  the  film  plane.  This  requires  accurate  calibration 
of  each  camera  and  accurate  determination  of  the  location  and  camera  orientation  in 
the  laboratory.  At  the  NBDL  an  optical  bench  with  a  plane  of  photo  targets  and  an 
Azimuth/Elevation  pedestal  for  the  canera  is  dedicated  to  the  determination  of 
individual  camera  calibration  parameters.  It  is  important  to  realize  that  there 
may  be  significant  differences  between  cameras. 

Low  frequency  errors  due  to  camera  orientation,  position  or  calibration  constants 
result  in  sharp  discontinuities  in  the  least  square  photo  solutions  when  the  target 
configuration  changes  due  to  drop-out  or  additions  of  targets.  These  discontinuities 
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are  accentuated  in  areas  where  the  solution  is  sensitive  to  errors.  In  order  to 
hold  these  discontinuities  to  a  minimum,  the  errors  in  determination  of  camera  calibration 
constants,  and  camera  location  and  orientation  should  be  minimized.  In  addition,  a 
computation  procedure  for  smoothing  the  discontinuities  may  be  required.  An  algorithm 
that  is  based  on  a  weighted  average  of  past  and  current  target  residuals  has  been 
developed  at  the  NBDL  to  remove  these  discontinuities. 

COMPARISON  OF  PHOTO  AND  ACCELEROMETER  DERIVED  ANGULAR  VELOCITY 

The  accelerometer  derived  variables  at  the  NBDL  are  derived  from  a  six  accelerometer 
configuration  with  sensitive  axes  orientation  and  accelerometer  location  as  shown  in 
Figure  2  (6).  The  six  accelerometer  system  has  been  much  maligned  in  favor  of  redundant 
nine  accelerometer  configurations.  Figures  3  through  6  present  angular  velocity 
around  the  anatomical  Y  axis  for  an  increasing  sequence  of  peak  G  levels  (5G  to  10G) 
for  -Gx  runs  derived  from  a  six  accelerometer  configuration  arranged  as  shown  in 
Figure  2. 

The  problem  is  three  dimensional,  but  only  the  significant  component  of  angular 
velocity  for  these  -Gx  runs  is  shown.  Figures  7  through  18  are  a  continuation  of 
this  sequence  for  13G  and  14G  peak  sled  acceleration.  A  three  triad  accelerometer 
configuration  was  used  for  these  runs.  The  three  triads  were  formed  by  adding 
accelerometers  where  required  to  make  each  triad  similar  to  the  triad  shown  in 
Figure  2.  Three  dimensional  components  of  angular  velocity  were  determined  both 
from  a  least  square  algorithm  for  the  three  triad  nine  accelerometers  and  the 
six  accelerometer  configuration  shown  in  Figure  2. 
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Figures  7  and  8  compare  the  angular  velocity  component  around  the  anatomical 
Y  axis  derived  from  the  three  triad  (nine  accelerometers)  and  six  accelerometers 
respectively  for  the  13G  run.  This  is  the  significant  component  of  angular  velocity 
for  the  -Gx  runs  and  as  seen  from  these  figures  very  little  difference  exists  between 
the  six  and  nine  accelerometer  solutions  and  both  agree  well  with  the  photo  derived 
angular  velocity.  A  similar  comparison  is  shown  for  the  angular  velocity  components 
around  the  anatomical  X  axis  in  Figures  9  and  10  and  the  anatomical  Z  axis  in 
Figures  11  and  12.  These  same  comparisons  are  shown  for  a  14G  run  in  Figures  13 
through  18.  The  solution  for  the  nine  accelerometers  is  better  than  that  for  the  six 
accelerometers  for  the  small  angular  velocity  components  around  the  anatomical 
X  and  Z  axes.  This  is  as  expected  since  there  is  a  statistical  advantage  with 
the  nine  accelerometers  and  the  error  in  either  case  is  significant  relative  to  the 
magnitude  of  the  angular  velocity  for  these  components  in  a  -Gx  run.  It  should  be 
noticed  that  the  profiles  for  the  six  and  the  nine  accelerometer  solutions  are  very 
similar  and  there  is  no  evidence  of  unstable  "blow  up"  in  either  case. 


SIMULATION  STUDIES 

A  program  has  been  developed  at  the  NBDL  to  evaluate  the  effect  of  low  frequency  accelerome 
errors  (sensitivity,  orientation  and  linearity)  on  the  errors  in  kinematic  variables 
derived  from  these  measurements.  This  program  uses  the  solution  of  the  kinanatic 
variables  to  generate  realistic  errors  in  each  of  the  accelerometers.  The  solutions 
obtained  with  the  accelerometers  modified  by  the  errors  are  compared  with  those 
obtained  from  the  unmodified  accelerometer  measurements.  In  this  way,  the  sensitivity 


to  errors  of  accelerometer  derived  kinematic  variables  is  evaluated  consistent  with 
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the  type  of  acceleration  profile  being  studied  at  the  NBDL.  Investigation  of  this  sort 
for  (-X),  (+¥)  and  (-X,  +Y)  vector  direction  profiles  studied  at  the  NBDL  have  indicated 
that  for  errors  in  accelerometers  consistent  with  pre-s election  '  accelerometers 
and  the  careful  calibration  procedure  there  is  no  blowup  of  sc  ions  with  six 

accelerometers.  The  errors  in  kinematic  variables  in  the  time  iow  of  interest 

(250  milliseconds)  are  consistent  with  the  error  growth  due  to  ^rating  the 
errors  in  angular  acceleration  resulting  from  deriving  angular  acceleration  components 
from  linear  acceleration  measurements.  These  conclusions  are  reinforced  by  the  good 
agreement  between  photo  and  accelerometer  derived  angular  velocity  components  in  the 
actual  data  base.  The  simulation  studies  show  that  with  significantly  larger 
accelerometer  errors  unstable  behavior  with  six  accelerometers 

is  observed.  Under  these  conditions  the  error  grows  out  of  bounds  in  very  short 
times.  This  behavior  has  never  been  observed  in  any  of  the  human  runs  at  the  NBDL. 
Simulation  studies  comparing  three  triad  solutions  with  six  accelerometer  solutions 
have  indicated  that  when  the  errors  in  the  accelerometers  are  large  enough  to  cause 
unstable  behavior  of  the  solutions  in  the  time  window  of  interest  with  six  accelerometers 
the  solutions  with  the  three  triad  configuration,  although  stable,  is  of  poor  quality. 


COMBINATIONS  OF  PHOTO  AND  ACCELEROMETER  DATA 


As  mentioned  previously,  it  has  always  been  the  intention  to  combine  photo  and 
accelerometer  solutions  into  one  "optimum”  solution  that  would  provide  a  consistent 
set  of  variables  from  the  acceleration  to  displacement  levels  suitable  for 
modeling  and  mannikin  construction. 
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The  solution  from  photo  data  although  good  for  displacement  and  angular  orientation, 
requires  extensive  signal  conditioning  in  order  to  derive  angular  and  linear 
acceleration  variables.  In  addition,  there  are  often  gaps  in  the  photo  solutions 
due  to  inadequate  target  coverage.  In  view  of  this,  an  algorithm  has  been  devised  to 
use  the  photo  displacement  and  orientation  results  to  predict  errors  in  the 
sensitivity  and  orientation  of  each  accelerometer.  This  solution  uses  both  the  photo 
derived  results  and  the  apriori  accelerometer  error  components.  The  solution  is  then 
obtained  using  accelerometer  data  modified  as  predicted  with  the  above  algorithms.  A 
complete  description  of  this  algorithm  will  be  presented  in  a  future  paper. 

Figures  19  through  24  compare  the  linear  and  angular  orientation  -.omponents 

derived  from  photo  measurements,  a  six  accelerometer  configuration  and  the 

optimum  combination  of  photo  and  accelerometers.  These  results  are  presented 

for  a  -7Gx  accelerometer  profile.  The  (+),  (X)  and  (A)  plot  symbols  are  consistent 

with  the  photo,  accelerometer  and  combined  solution,  respectively.  As  can  be 

easily  observed,  the  combined  solution  is  indeed  very  close  to  the  solution 

derived  from  the  photo  data.  Angular  velocity  comparisons  are  made 

in  Figures  25  through  27  for  this  same  -Gx  run.  Changes  in  linear  and  angular 

acceleration  components  are  obtained  but  are  small  and  are  not  shown.  The 

method  appears  encouraging  and  the  authors  intend  to  more  fully  explore  its  us;s  and 

implementation  in  the  future. 
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CONCLUS IONS 


1.  With  pre-selection  of  accelerometers  and  careful  calibration  procedures  six 
accelerometer  configurations  can  provide  results  which  agree  well  with  photo- 
derived  results  at  the  velocity  level  in  the  time  window  of  interest  (250 
milliseconds)  for  the  acceleration  profiles  investigated  at  the  NBDL.  The  error 
growth  in  the  time  window  of  interest  is  consistent  with  deriving  angular 
acceleration  from  linear  acceleration  measurements  and  integrating  the  low 
frequency  errors  in  the  accelerometers. 

2.  Given  the  same  accuracy  in  accelerometers,  there  is  a  statistical  advantage  of 
a  three  triad  (nine  accelerometer)  configuration  over  a  3-2-1  six  accelerometer 
configuration  and  the  nine  accelerometer  solution  is  marginally  stable.  However, 
if  the  errors  in  the  accelerometers  are  such  that  the  six  accelerometer  solution 
exhibits  unstable  behavior  in  the  time  window  of  interest  for  acceleration  profiles 
run  at  the  NBDL,  then  the  three  triad  solution  is  of  poor  quality. 

3.  An  algorithm  has  been  developed  which  combines  photo  and  accelerometer  data 
into  one  consistent  data  base.  This  algorithm  uses  a  model  for  accelerometer 
errors  and  can  use  non-contiguous  photo  derived  variables  to  predict  the  parameters 
in  these  models.  The  kinematic  variables  are  then  derived  from  the  accelerometer 
measursnents  and  the  modified  parameters  associated  with  each  accelerometer. 


— 
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Figure  1  T-Plate  Plots  for  Editing  Photo  Data 

Figure  2  Accelerometer  Configurations 

Figure  3  Angular  Velocity  Around  Anatomical  Y  Axis 
(-5Gx,  6  Accelerometers) 

Figure  4  Angular  Velocity  Around  Anatomical  Y  Axis 
(-7Gx,  6  Accelerometers) 

Figure  3  Angular  Velocity  Around  Anatomical  Y  Axis 
(-9Gx,  6  Accelerometers) 

Figure  6  Angular  Velocity  Around  Anatomical  Y  Axis 
(-1 OGx,  6  Accelerometers) 

Figure  7  Angular  Velocity  Around  Anatomical  Y  Axis 
(-13Gxf  9  Accelerometers) 

Figure  8  Angular  Velocity  Around  Anatomical  Y  Axis 
(-13Gx,  6  Accelerometers) 

Figure  9  Angular  Velocity  Around  Anatomical  X  Axis 
(-13Gx,  9  Accelerometers) 

Figure  10  Angular  Velocity  Around  Anatomical  X  Axis 
(-!3Gx,  6  Accelerometers) 

Figure  11  Angular  Velocity  Around  Anatomical  Z  Axis 
(-13Gx,  9  Accelerometers) 

Figure  12  Angular  Velocity  Around  Anatomical  Z  Axis 
(-13Gx,  6  Accelerometers) 

Figure  13  Angular  Velocity  Around  Anatomical  Y  Axis 
(-l4Gx,  9  Accelerometers) 

Figure  14  Angular  Velocity  Around  Anatomical  Y  Axis 
(-l4Gx,  6  Accelerometers) 

Figure  15  Angular  Velocity  Around  Anatomical  X  Axis 
(-l4Gx,  9  Accelerometers) 

Figure  16  Angular  Velocity  Around  Anatomical  X  Axis 
(-l4Gx,  6  Accelerometers) 


Figure  17  Angular  Velocity  Around  Anatomical  Z  Axis 
(-l4Gx,  9  Accelerometers) 


Figure  18 

Figure  19 

Figure  20 
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Figure  22 
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Figure  24 
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Figure  27 


Angular  Velocity  Around  Anatomical  Z  Axis 
(-14Gx,  6  Accelerometers) 

Comparison  of  Displacement  of  Head  Anatomical  Origin 
Relative  to  the  Sled  (Laboratory  X  Direction) 

Comparison  of  Displacement  of  Head  Anatomical  Origin 
Relative  to  the  Sled  (Laboratory  Y  Direction) 

Comparison  of  Displacement  of  Head  Anatomical  Origin 
Relative  to  the  Sled  (Laboratory  Z  Direction) 

Comparison  of  Ruler  Angle  Around  Head  Anatomical  X  Axis 

Comparison  of  Ruler  Angle  Around  Carried  Head  Anatomical  Y  Axis 

Comparison  of  Euler  Angle  Around  Carried  Head  Anatomical  Z  Axis 

Comparison  of  Angular  Velocity  Around  Head  Anatomical  X  Axis 

Comparison  of  Angular  Velocity  Around  Head  Anatomical  Y  Axis 

Comparison  of  Angular  Velocity  Around  Head  Anatomical  Z  Axis 
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'SECONDS  18  JUN  80 


RUN  NO.  SYM.  PLOTTED  SUBJECT  DELTA- 

LX3538  DAXSOA  =  CAD  H00108  .0020 

LX3538  DAXSOP  =  (XD  H00108  .0022 

LX3538  DAXSOS  =  (+D  H00108  .0020 


Sd313W  -  lN3N33VldSIQ 


RUN  NO.  SYM.  PLOTTED  SUBJECT  DELTA-T  ACC.  ONSET 
LX3538  DAYSOA  =  CAD  H00108  .0020  71.1  1561. 
LX3538  DAYSOP  =  C  X  3  H00108  .0022  71.1  1561. 
LX3538  DAYSOS  =  C+3  H00108  .0020  71.1  1561. 


TIME  -  SECONDS  23  MAY  80 


RUN  NO.  SYM.  PLOTTED  SUBJECT  DELTA-T  ACC 
LX3538  PHB02A  =  (A3  H00108  .0020  71. 
LX3538  PHB02P  =  (X)  H00108  .0022  71. 
LX3538  PHB02S  =  (+3  H001G8  .0020  71. 
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RUN  NO.  SYM.  PLOTTED  SUBJECT  DELTA-T  ACC.  ONSET 

LX3538  ■  PHC03A  =  (A)  H00108  .0020  71.1  1561. 

LX3538  PHC03P  =  (X)  H00108  .0022  71.1  1561. 

LX353S  PHC03S  =-  C  +>  —  HOftl  08-  .0020  71.1  1561. 
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RUN  NO.  SYM.  PLOTTED  SUBJECT  DELTA-T  ACC 
LX3538  RHAOXA  =  CAD  H00108  .0020  71. 
LX3538  RHAOXP  =  CX)  H00108  .0022  71. 
LX3538  RHAOXS  =  (+)  H00108  .0020  71. 
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RUN  NO.  SYM.  PLOTTED  SUBJECT  DELTA- 

LX3538  RHCOXA  =  (A )  H00108  .0020 

LX3538  RHCOXP  =  CX)  H00108  .0022 

LX3538  RHCOXS  =  (+)  H00108  .0020 


